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Abstract. We give a local central limit theorem for simple random walks 
on Z , including Gaussian error estimates. The detailed proof combines stan- 
dard large deviation techniques with Cramer-Edgeworth expansions for lattice 
distributions. 



A simple random walk on Z d is a sequence of independent identically distributed 
random variables (steps) on Z rf . We here suppose that the steplength is bounded, 
that is, the measure describing a single step has bounded support, and our goal is 
to give local approximations for the sum of the first n steps of the walk. Of course 
the Central Limit Theorem applies, and there also exist local CLTs, see for example 
Bhattacharya/Rao, |Q. 

In our work on local approximations for polymer measures (see the joint paper 
with E.Bolthausen |!|), we considered weakly self-avoiding walks as perturbations 
of simple walks. In order to obtain good bounds for the perturbations we thus 
needed good bounds for the simple walks as well. In particular, we needed Gaussian 
tail estimates, which were not provided by [El. Therefore we proved a version of a 
local CLT, which was tailored for our application there. 

It turns out that such a LCLT is useful not just in the treatment of polymer 
measures, but in all cases where similar perturbation methods are used and hence 
good local bounds for simple random walks are needed. 

Therefore we present a local central limit theorem with Gaussian tail estimates 
for walks with bounded steplength. This restriction may exclude some possible 
applications, but the setting is sufficient for many purposes and it allows a pre- 
cise control on the error terms. A major advantage over JIJ is that we make the 
dependence of the constants explicit. 

In order to formulate the theorem, we need some definitions and notations: Let 
G be the single step distribution of a random walk on Z d . The walk has bounded 
steplength i if G(x) — for all \x\ > i. The walk is called maximal, if the support 
of G is not contained in any affine hyperplane of 7L d . The walk is called aperiodic, 
if the greatest common divisor of all n € N with G* n (0) > is one. 

In the whole paper, tpr denotes the density of the centered normal distribution 
on M. d with covariance matrix T, that is, 

ipr{x) = (2 7 r)-' i / 2 (detr)- 1 / 2 exp[-x • {T' 1 ■ x)/2}. 

Furthermore we abbreviate if( n -id d ) by </?r ( , that is 

¥>„(«) = (27r V )- d / 2 exp[-x 2 /(2? ? )]. 

Note that for x £ M. d we abbreviate \x\ k by x k . 
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Local Central Limit Theorem. Let G be the single step distribution of a maxi- 
mal and aperiodic random walk on 1 d with mean zero and bounded steplength £ G N. 
Let r denote the covariance matrix of G. Then there exist polynomials P3 and Pq 
of degree three and six, respectively, and a positive constant C such that for all x 
in1 d , a £ (0, 1/2) and for all natural n, 

\G* n (x)-[l + n- 1 / 2 P 3 (x/V^)+n- 1 P 6 (x/V^)]^Pnr(x)\ < Cn^^^^x). 

The coefficients of the polynomials depend (polynomially) on T~ l and on the mo- 
ments of G up to order four, whereas C can be chosen independently of the specific 
law of G, depending only on d, i, a and on a lower bound for the smallest eigenvalue 
ofT. 

Before proving the theorem, we make some comments on the assumptions. The 
fact that the measure has bounded support is essential for the proof and can not 
be given up easily. The other assumptions are only technical: By shifting the 
measure we can arrange zero expectation and by restricting to a subspace we achieve 
maximality. For periodic measures we have to adapt the inversion formula for 
lattice measures ([To]). The effect of the period is that for fixed n we have non zero 
measure of G* n only on a sublattice of Z d , which leads to an additional factor in 
the approximating term but leaves the rest of the proof unchanged. 

The proof of the LCLT is split in two parts. The first part contains the large 
deviation arguments. In the second part we prove the Cramer-Edgeworth expan- 
sions, which are used in the first part to approximate the measures obtained by 
tilting the original distribution. 

Proof of the LCLT. Let 

Z(t) d = exp[i • x]G(x) 

and 1(0 = sup {t ■ £ - log Z(t)}. 

Standard large deviation theory (see for example (3j) yields a large deviation 
principle with entropy function / for the laws of G* n (x/n). Let Sg denote the 
convex closure of the set of points with nonzero G measure. Then / is convex on 
K d and even strictly convex on mt Sg, that is, the interior of Sg- Outside Sg, I 
equals +00. 

The function t 1— > D log Z(t) is an analytic diffeomorphism from M. d onto int Sg 
(for a proof see 13], page 261). Therefore, for any £ S int Sg, there exists a unique 
t 6 6 R d with DlogZ(i e ) = £. Clearly DlogZ(O) = and D 2 logZ(0) = V. For 
£ € int Sg, we have /(£) = t$ ■ £ - log Z(t^). Now for t <G R d denote by G t the tilted 
measure 

def G(x) exp[£ • x] 



Gt(x) 



Z(t) 



Using this, we see that for ^ d = x/n £ int Sg, we can write 
(1) G* n (x) = exp[-n/(£)] • Gl"(x). 

Case \x\ < n( 5 - Q )/ 9 : 

Since G is centered and maximal, the boundary of Sg is bounded away from 
zero. The covariance matrix of Gt ( depends analytically on £ with To = T. Let 
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7 be the smallest eigenvalue of T. We can find a constant 6 > such that for all 
£ G M d with |£| < 6 we have 

• £ £ int Sg and 

• the smallest eigenvalue of is greater or equal 7/2. 

We will now argue that S can be chosen depending only on t and 7, but not on 
the specific measure G. There is only a finite number of possible supports for a 
maximal and centered measure H with steplength f on Z d , and the convex hulls 
of these supports all contain zero in their interior. Therefore we find a constant 6' 
only depending on £, such that for all |£| < 5' we have £ S int for all such H. 

On the other hand we consider the set TL of all measures on Z d with steplength 
£, whose covariance matrices have smallest eigenvalue greater or equal 7. TL is 
obviously compact, and the measures in TL are all maximal (the covariance matrix 
of a non maximal measure has zero as an eigenvalue). The function that maps 
(H, £) to the smallest eigenvalue of the covariance matrix of H t( is continuous in 
both arguments. Hence we find a constant 6 < 6' such that the smallest eigenvalue 
of the covariance matrix of H ti is greater or equal 7/2 for all H £ TL and |£| < S. 
This shows the claimed independence. 

Now we come back to the proof of the LCLT. We know |£| = \x/n\ < ?i - ( 4 + Q )/ 9 
for a fixed constant a € (0,1/2). Thus for almost all n we have |£| < 8. ft is 
sufficient to prove the estimate for these n, since there is only a finite number 
of n and x with Sn < \x\ < n' 5_a ^ 9 . We can cover these cases by choosing C 
large enough. Moreover, using the compactness of TL again, we see that we find C 
depending on d, £, 7 and a only. 

So let (eS G with |£| < S A n~( 4+Q )/ 9 . To estimate the first factor in (Q), we 
use Taylor expansion for / at zero. 

Remember /(£) = if • £ — log Z(t^). We have DI(£) — and a simple computa- 
tion yields D 2 /(£) = T^ 1 . In general £) fc 7(C) (£,..., is a homogeneous polynomial 
of degree k, whose coefficients depend (polynomially) on T^ 1 and on the first to 
A:th moments of G^ . 

Now let TW denote a homogeneous polynomial of degree k. Using 1(0) — and 
DI(0) = 0, we obtain 



To deal with the error term R5 observe that, since the steplength of the tilted 
measures is also bounded by t, we can find upper bounds (depending only on d and 
€) for their moments. Furthermore we have • £| < |£|/7- Thus we know 



exp[-n/(0] = exp[-n£ • (T' 1 • 0/2] x exp[-,^ 3 ) (£) - nT^(i) - nO((£) 5 ))] 
= exp[-n£ • (r- 1 • 0/2] x [l - nT^(0 - nT^(0 + n 3 T< 6 >(0 



i(0 = £ ■ (r- 1 ■ 0/2 + (0 + t^ 4 \o + r 5 (i, 0, 0- 



\R 5 (I,0,0\ <const(c^,7)|£| 5 , 



and hence 



+ 0(n|£| 5 ) +Q(n 2 |£| 7 ) +Q(n 3 |£| 9 )exp[OH£| 3 )]] 



=0(n- 1 -°), since |f |<n-( 4 + = >/ 9 



(2) 



exp[-x- (T" 1 ■x)/(2n)\ 



n 



i (T (4) +T W)(x/y/n) + 0(- 



n 



-l-a 
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The coefficients of the poiynomials TW are depending polynomially on and 
the moments of G up to order four, and |0(n -1 ~")| < const (d, £, 7, a) nT 1 ^ . 

As a next step we approximate the second factor in (|l|), G*™(x). This is done 
in the second part of the paper, using Cramer-Edgeworth expansions for lattice 
distributions. Applying the corollary on page |§] on the measures Gt 5 , we obtain 

|G™(x)-(2vm)- d / 2 (det T^ 1 / 2 [l + n _1 L(0] | < C n -^ d+3 ^ 2 , 

where L(£) depends polynomially on the moments of Gt 5 up to order four and on 
the matrix T^ 1 . Since 7/2 is a lower bound for the smallest eigenvalue of T^, the 
constant C can be chosen depending only on d, £ and 7. 

We know that T^, T^ 1 and the other moments of Gt ( are analytic functions of 
£. Taylor expansion around £ = thus yields 

G%(x) = (27rn)- d / 2 (detr)- 1 / 2 

x [1 + T«(0 + T< 2 )(0 + 0(£ 3 ) + n- 1 ^") + 0(0] + 0(n" 3 / 2 )] 

(3) = (27rn)- d / 2 (detr)- 1 / 2 

x [1 + n-^ 2 T^-\x/y/n) + nT^T® + T^){x/y/n) + Ofr" 1 -")] , 

with the same dependencies as in (||). Inserting (@) and (||) in (Q) yields 

G* n (a:) = <p nr (x) x [l - r^^P^O/Vn) + n' 1 P {6) {x / Jn) + Ofa -1-0 ')] , 

where P^ denotes a polynomial of order k. Here P^ contains only odd order 
terms and P^ only even order terms. 

We still have to estimate (f n r- Since the coefficients of the covariance matrix T 
are bounded by £ 2 , we can bound the maximal eigenvalue of T by d£ 2 . Therefore 

ip n r(x) = {2irn)- d/2 {detry 1/2 exp[-{x ■ T- 1 x)/(2n)] 

< (2vm7)- d / 2 cxp[-a; 2 /(2^ 2 )] 

< (df h) d ' 2 W {x) 

< (2d£ 2 / 1 ) d / 2 Vn2M 2(x). 

Choosing C large enough yields the claim for \x\ < n^ a ^ 9 . 
Case \x\ > n( 5 - Q )/ 9 : 

For 'big' x we estimate G* n (x) and [l + n~ x l 2 Pz(x j yfn) + nT 1 Pq(x / ^fn)\ ip n r(x) 
separately. For fixed natural k we have \x\ k exp[— x 2 ] < const (fc) exp[—x 2 /%/2] for 
all x G Z d . Therefore |l + n~ l / 2 P^{x/ \/n) + n _1 Pq(x / \fn)\ip n v{x) is bounded by 
const(d, £, 7) f^/2 n ,i(2 (x) ■ Now we use the fact that there exists a natural number k 
such that n 1+a < (x/y/ri) k (k depending on a). Thus we find C = const(d,£,-f,a) 
such that for all n 

(4) |l + n- 1/2 P 3 (x/yf7i) + n- 1 PQ{x/y n R)\(p nT {x) < Cn~ l ~ a ip2ndp{x)- 

Now consider G* n (x). If £ ^ Sg, G* n (x) equals zero. If £ G int Sg, we use Taylor 
expansion of to obtain 

(5) 1(0= f (i-s)D 2 i(sO(Z,Ods= / 1 (i-,s)(e.r; ? 1 .£)d s >e 2 /(2^ 2 ). 

Jo Jo 
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In the last inequality we used that the maximal eigenvalue of T s ^ is majorized by 
d£ 2 for all £ and s. Inserting (§) into (§ yields G*"(jc) < exp[-x 2 / {2d£ 2 n)]. Now 
the same arguments that lead to (|j) yield the desired estimate for £ G int . 

Finally we consider £ £ dSc- Approximating £ with a sequence (£fc)fceN of points 
in int Sg, the corresponding sequence of measures Gt (k converges to a measure G'^ 
with support on a subset of the support of G, and — > /(£)■ Using (|l]) and (||) 
we therefore obtain 

G* n {x) = exp[-nl(0] • G^ n (x) < exp[-a; 2 /(2^ 2 n)]. 

The rest of the argument proceeds as before. □ 

CRAMER-EDGEWORTH EXPANSIONS FOR LATTICE DISTRIBUTIONS 

This second part yields the bounds we use for the tilted measures in the LCLT. 
The proof follows the arguments in Bhattacharya/Rao, Q|. The statement here 
differs from their result in two points: On one hand, we only want to approximate 
the walk at time n up to order n~( d+3 ^ 2 , and we suppose bounded range of the 
measure. Thus we only look at a special case of their result. On the other hand, we 
need more precise control of the error terms. This means that most of the required 
lemmas from Q would have to be adjusted. Rather than doing that, we prove 
everything we need directly. 

The object of interest here is H, the single step distribution of an aperiodic and 
maximal random walk on Z d with bounded steplength £ E N. Let E denote the 
expectation and Y the covariance matrix of H. The aim is to give a good local 
bound for H* n (nE), which we use in the proof of the LCLT. 

First we introduce some notations and basic estimates (for more details see 
chapter 6). We denote by H the Fourier transform of H , that is 

H(t)= e lt ' x H(x). 

Since H has bounded steplength, all moments of H exist. We denote by [i v the i^th 
moment of H, where v is a d-dimensional integral vector, that is 

M „ d ^ f x u H(x), and we have i H /i„ = (D v H)(0). 
The so called cumulants \v are given by 

and can be expressed in terms of moments by equating coefficients of t" on both 
sides of the formal identity 

\u\>l s=l \u\>l 

For r € N we define the polynomial Xr( z ) in z G C d by 

Xr{z) = ^ XvZ u . 

\u\=T 



6 



CHRISTINE RITZMANN 



We have xi(z) = E ■ z and Xii z ) = z ■ (T ■ z). For r S N we can bound 

Xr(z) 



|j/|=r 



Nfe+ll 

; k 



E 

,^ i )i,...,i„< fc )i>i 



^(i) . . . fj, v (_k) —j: — 



<M^ r E^(E i) fe <N r r[£(e d -i)r 

fe=l ll/l>l 



(7) <r(*f|z|r, 



where if ^ i[e d — 1). Thus for |i| < 1/K the Taylor series X^^Li ^fr^ is convergent 



00 Xr(tt) 



that is, 
(8) 



r=3 



Lemma. There exists a constant C such that for all x in Z d , a e (0, 1/2) and for 
all natural n, 



H* n {x)-n- d/2 



X3 (D) Xa{D) xf{D) 
6 y/n 24 n 72 n 



2 < T) > fx-nE 



(9) < C n-( d+3 )/ 2 , 

where the formal polynomials Xr(D) are given by 



Xr 



(D)f(x) ^ J2 ~Xu D"f{x) and x* r 2 (D)f(x) ^ ( Xr (D)f * Xr(D)f)(x) 



\y\=r 



The constant C depends only on the dimension d, the steplength t of H and on a 
lower bound for the smallest eigenvalue ofT. 

Proof. First we use the Fourier inversion formula for aperiodic lattice measures (see 
for example M (21.28)) to express H. We have 

H* n (x) = (2n)- d f exp[-it ■ x] H n (t) dt 



(10) 



-d 



X — t 

exp[— it ■ —= + n\ogH(— =)] dt. 



With the usual inversion formula for continuous measures we obtain 
(11) n- d ' 2 



_ X3 (D) X4(D) xf(D) ] ( x-nE ^ 



(2nVn-)- d I exp^.^^-^^lh+^ + ^M + ^M 



dt. 



where we used 



D"<pr(t) = (-it)"£P(i) = (-it)"exp[- 



t-(T-t)- 
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Now let £ be a small ball around with radius e > (precise conditions on 
e = e(d,£,j) will be determined later). Using ([ll]) and ( |ll| ) we can bound the left 
hand side of (||) by the sum of three integrals /, II and III with 

^ ri f r ~ , t sl r nE t-(T ■t) 1 
I = (27r^)- d / exp[nlogtf(^)] - exp\it • -= V — -L] 



and 

II = {2-Ks/n)~ d 
and 

III = {2ir^)- d 



[ — y/n-K , ,/mr ] d \ y/n£ 



1 



dt, 



Xajit) Xijit) xj(it) 
&Jn 24 n 72 n 



dt. 



exp[- 



t-(r-t)- 



i + 



6 Jn 



24 n 



72 n 



dt. 



We start with the essential term, the integral /. As a first condition we let 
e < 1/(2*0. Then we have K\t\/s/ri < 1/2 for all t S s/nE. Denote by 7 (a lower 
bound for) the smallest eigenvalue of T. We use (||) and (Q) to obtain 

t-(r-*)i 



\/n£ 



exp[- 



ex p[E'' 



Xr(H)- 



r=3 



6Vn 



24 n 



72 n 



< (2^)^ f exp[-^ 2 ] x [£ 



|Xr(**)| 



\/n£ 
1 

2 



r— 5 



r,s>3 
r+s>7 



^|Xr(**)||X.(**)l 



OO OO 

•£(£' 

s— 3 r— 3 



< (27rVn) 



-d 



exp[-^ 2 ] n-^Y.riKty + -n- 3 / 2 £ r S (i«) r+s 



\fn£ 



t— 5 



r,s>3 
r+s>7 



+ ]Tn- s/2 (2!rt) 3 7 S ! 

< n-( d + 3 ^ 2 (2^)- d f exp[- J< 2 ] const(d, £) [t 5 + i 7 + i 9 exp[ £ (2*:) 3 )t 2 ]] 

< n -(d+ 3 )/2 ^-d /" const ( d) £ )7 ) exp[-^ 2 ]cft, 

if we choose e = e(d,£, r f) small enough. Therefore we have 

I < const(d,£,7) n- (d+3)/2 . 

Now we come to the error term II. Since \H\ = 1 in zero and \H\ < 1 away 
from zero, there exists a constant K' < 1 with |if(t)| < K' on [— 7r, 7rp \ £. As 
in the proof of the LCLT we use the fact, that the set of all measures on Z d with 
steplength I such that the smallest eigenvalue of the covariance matrix is greater 
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or equal 7, is compact. Since the mapping (H,t) i— > H(t) is continuous in both 
arguments, the constant K' can be chosen depending on d, £ and 7 only. Therefore 
we have 

II = {2ir)- d f \H n (t) \ dt < K' n < const(d,£,7) n~( d+3) / 2 . 

J[-1T,Tt\ d \£ 

III is not much more complicated. Using (Q) to majorize the \r, we can bound 



III = (2ir^/n)- d / exp 

JR d \^/n£ 



t-(r-t) , 1+X a(it) , X4(«) , X§(«) 



6Vri 24n 72n 

2 



< (2ny/n)~ d / const(d,£)exp[-^-] 

JR d \^/n£ 3 
2 

< const(d^, 7 )n" d / 2 exp[-^^] < const(d^, 7 ) n - {d+ ^l 2 . 

The combination of these three estimates completes the proof of the Lemma. □ 

For the proof of the Local CLT we only need an approximation of H* n (x) at the 
point x = nE. In this case the approximating term 

X 3 (£>) , xdP) , X* 3 2 (D)^ ix- nE^ 



n- d ' 2 



1 



6y« 24 n 72 n 
simplifies to 

n^ 2 [1 + 7i- 1 ( X4 (^)/24 + xf (D)/72)] ^p(O), 

since the odd derivatives of a centered normal density vanish at zero. Therefore we 
obtain the following corollary as a special case of the lemma. 

Corollary. There exist constants C and L such that for all natural n, 

\H* n {nE) - (2 7 rn)-' i / 2 (detr)- 1 / 2 [l + n^L] | < C n-^ d+3 ^ 2 . 

Here L is a constant depending polynomially on the moments of H up to order four 
(coming from the x terms) and on the matrix (coming from the derivatives of 
tpr )■ The constant C depends only on the dimension d, the steplength t of H and 
on a lower bound for the smallest eigenvalue ofT. 

References 

[1] Bhattacharya, R. N.; Rao, R. R.: Normal Approximation and Asympotic Expansions. 

Reprint, Malabar 1986 (Original 1976). 
[2] Bolthausen, E.; Ritzmann, C: A Central Limit Theorem for Convolution Equations 

and Weakly Self- Avoiding Walks. Preprint, Zurich 2002. 
[3] Ellis, R. S.: Entropy, Large Deviations and Statistical Mechanics. New York 1985. 



